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AND MULTIPLIER IDEALS 
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Abstract. The Bernstein-Sato polynomial (or global fe- function) is an important invari- 
ant in singularity theory, which can be computed using symbolic methods in the theory 
of Z?-modules. After providing a survey of known algorithms for computing the global 
6-function, we develop a new method to compute the local &-function for a single polyno- 
mial. We then develop algorithms that compute generalized Bernstein-Sato polynomials 
of Budur-Mustat;a-Saito and Shibuta for an arbitrary polynomial ideal. These lead to 
computations of log canonical thresholds, jumping coefficients, and multiplier ideals. Our 
algorithm for multiplier ideals simplifies that of Shibuta and shares a common subroutine 
with our local 6-function algorithm. The algorithms we present have been implemented 
in the Z?- modules package of the computer algebra system Macaulay2. 



1. Introduction 

The multiplier ideals of an algebraic variety carry essential information about its singu- 
larities and have proven themselves a powerful tool in algebraic geometry. However, they 
are notoriously difficult to compute; nice descriptions are known only for very special fam- 
ilies of varieties, such as monomial ideals and hyperplane arrangements [T9 | |35 | |27] . To 
briefly recall the definition of this invariant, let X = C" with coordinates 
For an ideal (/) = (/i, • • • , /r) ^ C[a3] and a nonnegative rational number c, the multiplier 
ideal of f with coefficient c is 



(El/. 



2V 



is locally integrable 



It follows from this definition that J{f) D J{f) ioi c < d and J{f) = C[x] is 
trivial. The (global) jumping coefficients of / are a discrete sequence of rational numbers 
= with = ^0 < ^1 < ^2 < ■ ■ ■ satisfying the property that Ji^f^) is constant 

exactly for c G In particular, the log canonical threshold of / is ^i, denoted by 

lct(/). This is the least rational number c for which is nontrivial. The multiplier 

ideal J{f'^) measures the singularities of the variety of / in X; smaller multiplier ideals 
(and lower log canonical threshold) correspond to worse singularities. For an equivalent 
algebro-geometric definition and an introduction to this invariant, we refer the reader to 



In this paper we develop an algorithm for computing multiplier ideals and jumping coef- 
ficients by way of an even finer invariant, Bernstein-Sato polynomials, or ^-functions. The 
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results of Budur et al. [6] provide other applications for our Bernstein-Sato algorithms, in- 
cluding multiplier ideal membership tests, an algorithm to compute jumping coefficients, 
and a test to determine if a complete intersection has at most rational singularities. 

The first 6-function we consider, the global Bernstein-Sato polynomial of a hypersurface, 
was introduced independently by Bernstein [1] and Sato [29]. This univariate polynomial 
plays a central role in the theory of D-modules (or algebraic analysis), which was founded 
by, amongst others, Kashiwara [TT] and Malgrange [TTj. Moreover, the jumping coeffi- 
cients of / that lie in the interval (0, 1] are roots of its global Bernstein-Sato polynomial 
[?]; however, this 6-function contains more information. Its roots need not be jumping 
coefficients, even if they are between and 1 (see Example 16. ip . 

The Bernstein-Sato polynomial was recently generalized by Budur et al. |6] to arbi- 
trary varieties. The maximal root of this generalized Bernstein-Sato polynomial provides 
a multiplier ideal membership test. Shibuta defined another generalization to compute ex- 
plicit generating sets for multiplier ideals [52]. Our multiplier ideal algorithm employs the 
6-functions os Shibuta, which we call the m- generalized Bernstein-Sato polynomial. How- 
ever, it circumvents primary decomposition and one elimination step through a syzygetic 
technique (see Algorithms 14.51 and I3.2p . The correctness of our results relies heavily on 
the use of l^-filtrations, as developed by Kashiwara and Malgrange [121 HH] • 

D-module computations are made possible by Grobner bases techniques in the Weyl 
algebra. The computation of the Bernstein-Sato polynomial was pioneered by Oaku in 
[21]. His algorithm was one of the first algorithms in algebraic analysis, many of which 
are outlined in the book by Saito et al. [21]. The computation of the local Bernstein- Sato 
polynomial was first addressed in the early work of Oaku [21], as well as the recent work 
of Nakayama [20], Nishiyama and Noro [21], and Schulze [301 EI]- Bahloul and Oaku [2] 
address the computation of local Bernstein-Sato ideals that generalize Bernstein-Sato 
polynomials. In this article we provide our version of the local algorithm for Bernstein- 
Sato polynomials, part of which is vital to our approach to computation of multiplier 
ideals. 

There are several implementations of algorithms for global and local 6-functions in 
kan/sml [33], Risa/Asir [23], and Singular [9]. One can find a comparison of performance 
in [15] . All of the algorithms in this article have been implemented and can be found in 
the D-modules package [TB] of the computer algebra system Macaulay2 [S] . 

The first author was partially supported by NSF Grants DMS 0555319 and DMS 090112; 
the second author is partially supported by the NSF Grant DMS 0914802. 



Outline. Section |2] surveys the known approaches for computing the global Bernstein- 
Sato polynomial, highlighting an algorithm of Noro [22]- In Section [3l we present an 
algorithm for computing the local Bernstein-Sato polynomial. Algorithms for the gen- 
eralized Bernstein-Sato polynomial for an arbitrary variety, as introduced by Budur et 
al. [nj , are discussed in Section HI along with their applications. Based on the methods 
of Section [31 Section considers the m-generalized Bernstein-Sato polynomial of Shibuta 
[32] and contains our algorithms for multiplier ideals. 
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2. Global Bernstein-Sato polynomials 

Let K he a field of characteristic zero, and set X = and Y = X^K witli coordinates 
{x) and respectively. We consider the n-th Weyl algebra Dx = K{x,d) with 

generators Xi, . . . ,Xn and d^^, . . . , d^^, as well as Dy = K{x, dx, t, dt), the Weyl algebra 
on Y. Define an action of Dy on Nf := K[x][f~^, s]f^ as follows: Xi and act naturally 
for i = 1, . . . ,n, and 

t-hix,s)r = h{x,s + l)fr and dt ■ h{x, s)r = ~sh{x, s - l)r' f , 

where h G K[x][f~^, s]. 

Let a = —dtt. For a polynomial / G -ft'[ic], the global Bernstein-Sato polynomial of /, 
denoted 6/, is the monic polynomial h{s) G K[s\ of minimal degree satisfying the equation 

(2.1) b{a)r = pfr 

for some P G Dx{o'). 

There is an alternate definition for the global Bernstein-Sato polynomial in terms of 
y-filtrations. To provide this, we denote by V'Dy the ^-filtration of Dy along X, where 
V^Dy is Dx-generated by the set {f^d^ \ — v > m]. Let if : X Y defined by 
if{x) = {x, f{x)) be the graph of /. The D-module direct image of K[x] along if is the 
module 

Mf := {if)+K[x] = K[x] ®K K{dt) 
with actions of a vector field on X and t, 

^(p ® 5^ = ® - ® and t- {p^d'^) = fp0d'^ -up®d'^-\ 

providing a Dy-module structure. Notice that there is a canonical embed ding of Mf into 
Nf, where s is identified with —dtt. 

With (5 = 1 ® 1 G Mf, the global Bernstein-Sato polynomial 6/ is equal to the minimal 
polynomial of the action of a on the module {V^ Dy)6 / {V^ Dy)6 . We now survey three 
ways of computing this 6-function. 

2.1. By way of an annihilator. The global Bernstein-Sato polynomial bf{s) is the 
minimal polynomial of a := —dtt modulo Ann^i^-jo-] + Dx[(j]f, where G Nf. By the 
next result, this annihilator can be computed from the left Dy-ideal 

If = {t-f,^^ + ^^^t,...,^^ + §^^t). 

Theorem 2.1. [281 Theorem 5.3.4] The ideal Ann£)[s] equals the image of If n D[a] 
under the substitution a t— )■ s. 

2.2. By way of an initial ideal. This method makes use w = (0, 1) G x M, the 
elimination weight vector for X in Y. 

Theorem 2.2. Let b{x,s) be nonzero in the polynomial ring K[x,s]. Then b{x,a) G 
{i'n.(^_^^y^)If)r\K[x,a] if and only if there exists Q G D[s] satisfying the functional equation 
Qjs+i _ ^(^2;, s)^^. In particular, 

{bf{a)) = in(_^,^)// n K[a]. 
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Proof. The action of t on Nf is multiplication by /, hence, the existence of the functional 
equation is equivalent to b{x, s) G // + V^Dy- The result now follows from Theorem 12. 
which identifies s with a. □ 

The following algorithm provides a more economical way to compute the global b- 
function using linear algebra. By establishing a nontrivial i^'-linear dependency between 
normal forms NFg'(s*) with respect to a Grobner basis G of in(_^ where < i < d 
and d is taken as small as possible, this algorithm bypasses elimination oi di, . . . , dn- This 
trick was used for the first time by Noro in [22], where a modular method to speed up 
6-function computations is provided as well. We include the following algorithm for the 
convenience of the reader as a similar syzygetic approach will be used in Algorithms 13.21 
14.51 and 15. 12] Note that the coefficients of the output are, in fact, rational, since the roots 
of a 6-function are rational [TT] . 

Algorithm 2.3. b = global B Function{f , P) 
Input: a polynomial f G 

Output: polynomial b G Q[s] is the Bernstein-Sato polynomial of f . 
G ^ Grobner basis o/in(_^ „,)Jy. 
d^O. 
repeat 
d ^ d+1 

until 3{co, ■ ■ ■ ,Cd) G Q"'^-'^ such that q = 1 and 

d 

^c,NFg(sO = 0. 

i=0 

return Y.i^QCiS\ 

This approach can be exploited in a more general setting to compute the intersection 
of a left ideal with a subring generated by one element as shown in [T]. 

2.3. By way of Briangon— Maisonobe. This approach, which is laid out it [5], com- 
putes the annihilator of f^ in an algebra of solvable type similar to, but different from, the 
Weyl algebra. This path has been explored by Castro- Jimenez and Ucha |36] and imple- 
mented in Singular pj with a performance analysis given by Levandovskyy and Morales 
in [12] and recent improvements outlined in [T]. 

3. Local Bernstein-Sato polynomials 

In this section, we provide an algorithm to compute the local Bernstein-Sato polynomial 
of / at a prime ideal of i^[a;], which is defined by replacing the use of Dx in (12. ip by its 
appropriate localization. Algorithms 13.11 and 13.21 use Theorem 12.21 to compute an ideal 
Eh C K[x] that describes the locus of points where the 6-function does not divide the 
given b G Q[s]. 

Algorithm 3.1. Ef, = exceptionalLocusB(/, 6) 
Input: a polynomial f G -ft'fa;], a polynomial b G Q[s]. 
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Output: Eb C K[x] such that V P G SpecK[x], 

G ^ generators o/in(_^ fl K[x, s], where s = —dtt. 
return exceptionalLocusCore {G,b). 

The following subroutine computes /'^[ajj-syzygies between the elements of the form s^g 
of s-degree at most deg6 and b itself. It returns the projection of the syzygies onto the 
component corresponding to b. 

Algorithm 3.2. Ef, = exceptionalLocusCore (/, 6) 
Input: G G K[x, s], a polynomial b G Q[s]. 
Output: Eb C K[x]. 

Gi i— {a Grobner basis of (G) w.r.t. a monomial order eliminating s}. 

d ^ deg b. 

G2 ^ {s'g \ g eGi, i + deg, ^ < 4. 
S ^ ker where 

d 

(t) : K[x\\^^\+^ ^K[xy 

i=0 

maps Ci, for 2 = 1,..., \G2\, to the elements of G2 and eic^l+i '^'^ ^■ 
return projection of S C K[x]^'^^^~^^ onto the last coordinate. 

The computation of syzygies in line S] and projection in line Oof Algorithm 13.21 may be 
combined within one efficient Grobner basis computation. 

of correctness of Algorithms rO] and UTB . The local Bernstein-Sato polynomial bf^p at 
P e Spec [a;] divides the given 6 G Q[s] if and only if 

Q7^+i = bf^pf, for some Q' eK[x]p® D[s] 
^ Qf'+^ = hbf % for some Q e D[s], heK[x]\P. 

For h G K[x], 

Qf^+^ = hbf, for some Q E D[s] 
<^ hb e m^-u),w)If n K[x, s] (by Theorem [2^ 
<^ h is the last coordinate of a syzygy 

in the module produced by line H] 

h e Eb. 

This proves that bf^p\b^Eb(/LP. □ 

Remark 3.3. [Particulars of Algorithm 13. Ij In order to compute generators of in(_^^)Jj, 
one may apply the homogenized Weyl algebra technique (for example, see Algo- 
rithm 1.2.5]). Then to compute generators of in(_^_^)Jj fl K[x\{t,dt) , eliminate dx and 
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apply the map ip defined as follows: for a (— w, w) -homogeneous h G K[x]{t,dt) with 
deg(_^,^) h = d, 

f^K if > 0, 



This is the most expensive step of the algorithm 



d^K ifc^<0. 



We are now prepared to compute the local Bernstein-Sato polynomial of / at a prime 
ideal P C fTfcc]. Its correctness follows from that of its subroutine, Algorithm 13. II 

Algorithm 3.4. h = localBFunction(/, P) 
Input: a polynomial f G , a prime ideal P C 
Output: b G Q[s], the local Bernstein-Sato polynomial of f at P. 
b^bf. {global b-function} 
for r G bj^{0) do 
while [s — r) \ b do 
b' ^ b/{s-r). 

i/exceptionalLocusB(/, 6') C P then 

break the while loop, 
else 

b^b'. 
end if 
end while 
end for 
return b. 

Remark 3.5. Algorithm 13. II can also be used to compute the stratification of Speci^fa;] 
according to local 6-function. Below are the key steps in this procedure. 

(1) Compute the global fe-function bf. 

(2) For all roots c G bj^{0) compute 

Ec^i = exceptionalLocusB(6//(s — cY"'^), 

where i > and is at most the multiplicity fic of the root c in bf. 

(3) The stratum of = n^g^--i(Q)(s — c)*% a divisor of bf, is 



V 



n \ u ^(^' 



ycgb-i(o), jc>o J 



This approach is similar to that in the recent work [21] of Nishiyama and Noro, which 
offers a more detailed treatment. 

4. Generalized Bernstein-Sato polynomials 

4.1. Definitions. For polynomials / = /i, ...,/,• e K[x], let = 111= i and Y = 
X with coordinates {x,t). Define an action of Dy = K{x,t, d^, dt) on Nf : = 
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K[x][f~'^, s]f^ as follows: Xi and 9a;., for z = 1, . . . , n, act naturally and 

tj ■ h(^X, Si, . . . , Sj, . . . , Sr)f h{x^ Si, . . . , Sj • + 1, . . . , Sr)fjf , 

dt. ■ h{x, Si, . . . , Sj, . . . , Sr)f^ = -Sjh{x, Si, . . . , - 1, . . . , s^)//V^ 

for j = 1, . . . , r and /i G -R'[a:;] [Z^"*^, s]. 

With cr = — (X]i=i ^ti^j)? the generalized Bernstein-Sato polynomial bf^g of / at (7 G 
ii'[a;] is the monic polynomial b G C[s] of the lowest degree for which there exist Pk G 
Dx{dtitj \ 1 < i, j < r) for k = 1, . . . ,r such that 

r 

(4.1) Ka)gr = Y,P,ghr- 

k=l 

Remark 4.1. When r = 1, the generalized Bernstein-Sato polynomial bf i = bf is the 
global Bernstein-Sato polynomial of / = /i discussed in Section [21 

There is again an equivalent definition of bf^g by way of the V^-filtration. To state this, 
let V'Dy denote the ^-filtration of Dy along X, where V^^Dy is Z)x-generated by the 
set {t^d^ I — |z/| > m}. The following statement may be taken as the definition of the 
y-filtration on 

Theorem 4.2. [6, Theorem 1] Fore G Q and sufficiently small e > 0, = V^~^''K[x] 

and V^K[x] = J{f^'). 

Consider the graph of /, the map if : X defined by if{x) = {x, fi{x), . . . , fr{x)). 
We denote the D-module direct image of K[x] along if hj 

(4.2) Mf := {tf)+K[x] = K[x] ®k K{dt). 

This module carries a Dy-module structure, where the action of a vector field on X and 
that of tj are given by 

r 

aP ® ) = ® d'i - 5^(e/.)p ® ^r'^ and t, ■ ip ® d'i) = f)p ® - ^jP ® 9^"' , 

i=l 

where = YVi=i^ti ^ = i'^i, ■ ■ ■ ^i^r) ^ N and Cj is the element of with j-th 
component equal to 1 and all others equal to 0. 
Further, Mf admits a l^-filtration with 

V^Mf = J2 {V'''+^''^K[x]) ® dl 

For a polynomial g G K[x] so that (7 (8> 1 G Mf, bf^g is equal to the monic minimal 
polynomial of the action of a on 

_ {V'Dy)ig®l) 
{VWy){g®iy 

Remark 4.3. There is a canonical embedding of Mf into Nf, where Si is identified with 
—dt^ti. In particular, for a natural number m, the image of {y'"^Dy){l ® 1) under this 
embedding is contained in {V^Dy){f)"^ C Nf. 
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4.2. Algorithms. To compute the generahzed Bernstein-Sato polynomial, we define the 
left Dy-ideal 

// = {U-f^\l<t<r) + (9,, + E:=i H^x. I 1 < J < n) 
that appears in the following multivariate analog of Theorem 12.11 Recall that a = 

Theorem 4.4. The ideal If is equal to Arniuy f^. Furthermore, the ideal Ann^j^j^] 
equals the image of If fl -DxM under the substitution cr i— s. 

We now provide two subroutines in our computations of Bernstein-Sato polynomials 
and multiplier ideals. The first finds the left side of a functional equation of the form (14. ip 
without an expensive elimination step. The second finds the homogenization of a Dy-ideal 
with respect to the weight vector {—w,w), where w = (0, 1) G x W determines an 
elimination term order for X in Y. 

Algorithm 4.5. b = linearAlgebraTrick((7, G) 
Input: generators G of an ideal I C Dy, 

a polynomial g G -ft" [a?], 

such that there is b E K[s] with b{a)g G /. 
Output: b, the monic polynomial of minimal degree such that b{a)g G /. 

B ^ {a Grobner basis of DyG}. 

repeat 

d<- d+l 

until 3(co, . . . , Crf) G K''-^^ such that Q = 1 and 

d 

J2c^^FB{(J'g)=0. 
1=0 

return Yl'i=o^i^^- 
Algorithm 4.6. G* = starldea^G, w) 

Input: generators G of an ideal J C Dy, 

a weight vector w G Z""'"''. 
Output: G* C gTi^_^^-^Dy = Dy, a set of generators of the ideal J* of {—w,w)- 

homogeneous elements of J. 

G^ generators G homogenized w.r.t. a weight {—w,w); G^ C Dy[h] with a homoge- 
nizing variable h of weight 1. 

B ^ {a Grobner basis of (G'*, hu — 1) G Dy[h,u] w.r.t. a monomial order eliminating 
{h,u}}. 

return B fl Dy . 

Below are two algorithms that are simplified versions of Shibuta's algorithms for the 
generalized Bernstein-Sato polynomial. In the first, we use a module -Dy[s], where the 
new variable s commutes with all variables in Dy. 
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Algorithm 4.7. bf^g = generalB(/, (7, Starldeal) 
Input: / = {/i, . . . , M C K[x], g E K[x]. 

Output: bf^g, the generalized Bernstein-Sato polynomial of f atg. 
G,^ {tj'-f^\j = l,...,r} U {d,^ + ^^.^^^dt^ \ z = l,...,n}. 

G2 ^ starIdeal(Gi, w) U {gfi | 1 < i < r) U {s — a} C -Dy[s], where w assigns weight 1 
to all dtj and to all dx^ ■ 
return linear AlgebraTrick(G2)- 

Algorithm 4.8. bf^g = generalB(/, gf, Initialldeal) 
Input: f = {/i, ...Jr}C K[x], g E K[x]. 

Output: bfg, the generalized Bernstein-Sato polynomial of f atg. 
Gi ^ {t, - /, I J = 1, . . . , r} U {dx^ + E;=i g^*, h = 1, . . . , n}. 
G2 ^ GiH Dy ■ g. 

G3 ^ generators of 'm.(^_^ ^^{(^2) , where w assigns weight 1 to all dt^ and to all d^^. 
return linearAlgebraTrick(G3). 

Their correctness follows from [521 Theorems 3.4 and 3.5]. 

Remark 4.9. According to the experiments in [T3] a modification of Algorithm 14.61 that 
uses elimination involving one less additional variable exhibits better performance. Our 
current implementation does not take advantage of this. 

4.3. Applications. The study of the generalized Bernstein-Sato polynomial in [6] yields 
several applications of our algorithms, which we mention here. Each has been implemented 
in Macaulay2. 

We begin with a result that shows that comparison with the roots of bf^g{s) provides a 
membership test for Jif^) for any positive rational number c. 

Proposition 4.10. P, Corollary 2] Let g E K[x\ and fix a positive rational number c. 
Then g E Jif^) if and only if c is strictly less than all roots of bfg{—s). 

When / defines a complete intersection, Algorithms 14.71 and 14.81 provide tests to deter- 
mine if Z has at most rational singularities. 

Theorem 4.11. [HI Theorem 4] Suppose that Z is a complete intersection of codimension 
r in Y defined by f = fi, . . . , fr- Then Z has at most rational singularities if and only if 
lct(/) = r and has multiplicity one as a root ofbf{—s). 

To compute a local version of the generalized Bernstein-Sato polynomial, we need the 
following analog of Theorem 12. 2[ 

Theorem 4.12. Let b{x,s) be a nonzero polynomial in K[x,s]. Then the polynomial 
b{x,a) E in(„^^)Jj nK[x,a] if and only if there exist Qk E D[s] s.t. '^^^iQkfkf'^ = 

b{x,s)r- 

Proof. This follows by the same argument as that of Theorem 12. 2[ □ 

Remark 4.13. In light of Theorem I4.12[ the strategy in Section [3] yields a computation 
of the local version of the generalized Bernstein-Sato polynomial. The only significant 
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difference comes from the lack of an analogue to the map ip of Remark 13.31 However, it 
is still possible to compute 'm.{^-u},w)If H K[x, a] by adjoining one more variable s to the 
algebra and s — a to the ideal and eliminating t and dt- In case of the hypersurface this 
is a more expensive strategy than the one described in Remark 13.31 

5. Multiplier ideals via m-generalized Bernstein-Sato polynomials 

For this section, we retain the notation of Section H] and discuss Shibuta's m-generalized 
Bernstein-Sato polynomials. These are defined using the ^-filtration of Dy along X, but 
they also possess an equational definition. In contrast to the generalized Bernstein-Sato 
polynomials of Section HI this generalization allows us to simultaneously consider families 
of polynomials -^'[aj], yielding a method to compute multiplier ideals. 

Definition 5.1. Let mJ"^ := {V^ Dy)6 / {V^ Dy)5 with 5 = l®l e Mf = K[x\ ®k 
K{dt) ■ Define the m-generalized Bernstein-Sato polynomial b^^J to be the monic minimal 
polynomial of the action of a := — (X]i=i ^uti) on 

Remark 5.2. Since is ^-filtered, the polynomial b^J^^ is nonzero and its roots are 
rational. 

Proposition 5.3. The m-generalized Bernstein-Sato polynomial b^^g is equal to the monic 
polynomial b{s) of minimal degree in K[s] such that there exist Pk G Dx{—dt-tj \ I < i, j < 
r) and hk G (/)™' such that in Nf there is an equality 

r 

(5.1) h{a)gr = Y,Pkhkr- 

k=i 

Proof. By the embedding in Remark 14.31 the existence of such an equation is equivalent 
to the existence of Qk G Dx{—dt^tj \ I < i, j < r) and /x(/c) G N'" with |/i(fc)| > m such 
that in Mf, b{a) ■ {g I) = ELi Qkt''^''^ ■ (1 ® 1). □ 

Remark 5.4. Since (y^DY)g ® 1 C Mj'* is a quotient of Mf g, the generalized Bernstein- 
Sato polynomial bf^g is a multiple of the m-generalized Bernstein-Sato polynomial b^f^^. 

When g is a unit, the equality bf^g = fe^^ holds, as seen easily by comparing (14. ip and 
(15. ip . However, this equality does not hold in general. 

Example 5.5. When n = 3 and / = XlLi have 

bf,-iis) = {s + l){s + ^) and b^^l^{s) = s + 1 . 

In particular, b^f\-^ strictly divides 

Proposition 15.31 translates into the following algorithm. 
Algorithm 5.6. b^p^ = generalB(/, gf, m) 



ALGORITHMS FOR BERNSTEIN-SATO POLYNOMIALS AND MULTIPLIER IDEALS 11 

Input: f = {/i, ...,fr}c K[x], g e K[ 

Output: b^pg, the m- generalized Bernstein- Sato polynomial ( as defined in Definition \5. 1\) . 

Gi^ {t,-/, |j = l,...,r} U {9.^ + ^;^^|gajz = l,...,n}. 

G2 ^ starIdeal(G'i, ) U {/" | a G N'', |a| = m}, where w assigns weight 1 to all dt^ 

and to all ■ 

Jul 

return linear AlgebraTr ick( (7, G2) • 

5.1. Jumping coefficients and the log canonical threshold. For the remainder of 
this article, set K = C. Our algorithms for multiplier ideals are motivated by the following 
result. 

Theorem 5.7. [321 Theorem 4.3] For g G K[x] and c < m + \ct{f), g e Jif) if and 
only if c is strictly less than every root of bpg{—s) . In other words, 

JiD = {9^ K[x] I bp^{-a) = ^ c < «}. 

Proof. By Theorem 14. 2^ Jif^) = V'^^'^K[x\ for all sufficiently small e > 0. Hence, g G 
Ji^f^) precisely when g ®1 & V^Mf for all a < c, or equivalently, 

(5.2) c < max{a \g®le VMf}. 

As in P (2.3.1)], the right side of (jSJ]) is equal to min{a | Gr^((\/ODy)(c/ O 1)) ^ 0} 
and strictly less than min{a | Gry((y™Dy)5) 7^ 0}. Thus by our choice of m, g E J^f^) 

exactly when c is strictly greater than min{a | Gry M^™^ 7^ 0}. The theorem now follows 
because Gr^ M^]^ ^ if and only if b^^ {-a) =0. □ 

Theorem 15.71 provides a second membership test for membership in J{f^)\ moreover, 
the following corollary provides a method for computing the log canonical threshold and 
jumping coefficients of / via the m-generalized Bernstein-Sato polynomial b^ = b^^. 

Corollary 5.8. For any positive integer m, the minimal root of bp\—s) is equal to the 
log-canonical threshold lct(/) of (/) C K[x]. Further, the jumping coefficients of (/) 
within the interval [lct(/), lct(/) + m) are all roots of bp\—s) . 

5.2. Computing multiplier ideals. Here we present an algorithm to compute multiplier 
ideals that simplifies the method of Shibuta [32] • In particular, significant improvement 
is achieved bypassing the primary decomposition computations required by Shibuta's 
method. 

For a positive integer m, define the K[x, a]-ideal 

Jf{m) = {l} + Dy{fr)nK[x,a], 

where Ij^ C Dy is the ideal of the (— ty, u;)-homogeneous elements of If. This ideal is 
closely related to the m-generalized Bernstein-Sato polynomials. 
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Lemma 5.9. For g G the m- generalized Bernstein-Sato polynomial b^p^ is equal to 

the monic polynomial b{s) G K[s] of minimal degree such that 

(5.3) {bia)) = iJf{m):g)nK[a]. 

Proof. By (jEI]), fej"^^ is the monic polynomial b{s) G K[s] of minimal degree such that 

bia)g elf + Dx{-dti, I 1 < 2, J < r) ■ (/)™. 
Since b{a)g is (— w, w)-homogeneous, we obtain fl5.3p . □ 

Theorem 5.10. [321 Theorem 4.4] Let Jf{m) = flLi 9i be a primary decomposition with 
Qi n K[a] = (cr + c{i)Y^^^ for some positive integer K{i). Then for c < lct(/) + m, 

Jin= n (^.-n^w). 

j-cij)>c 

Proof. We see from fl5.1l) that b^Pg{s) is the monic polynomial b{s) of minimal degree 
such that there exist some Pk G Dx{—dtitj | 1 < i,J < t) and hk G such that 

{K'^)9 - J2k Pkhk) G //• Equivalent^, 

b{a)ge{rf + DY-{fr)nK\x,a]. 
The theorem now follows from Lemma 15. 9[ □ 



The following is based on methodology used in the computation of the local fe-function 
and, in particular, employs Algorithm 13.21 Its correctness follows immediately from The- 
orem 15.101 and the results of Section [31 

Algorithm 5.11. J{f^) = multiplierldeal(/, c) 

Input: / = {/i, ...,/.} C K[x], c G Q. 

Output: J'{f^), the multiplier ideal of f with coefficient c. 

Gi ^ {t, - fj I J = 1, . . . , r} U {9., + E;=i z = 1, . . . , n}. 

m ^ [max{c — lct(/), 1}] . 

if c — lct(/) is integer and > 1 then 
m ^ m + 1 

end if 

G2 ^ starIdeal(Gi, ty) U {/" | a G N*", |a| = m} U {s — a} C Dy[s], where w assigns 
weight 1 to all dt and to all dx^ • 

B ^ { a Grobner basis of G2 w.r.t. an order eliminating {dx,t,dt} } fl K[x,s]. 
b ^ generalB(/, l,m). {The computation ofb^^ may make use of B.} 

b' <r- product of factors (s — c')"*-'^'^ ofb over all roots d ofb^j^^ such that —c' > c, where 
a{c') equals the multiplicity of the root c' . 
return exceptionalLocusCore(i?, b'). 

As noted in [32l Remark 4.6.ii], Jif") = (/) • J{f^~^) when c is at least equal to the 
analytic spread A(/) of (/). (The analytic spread of (/) is the least number of generators 
for an ideal / such that (/) is integral over /.) Hence, to find generators for any multiplier 
ideal of /, it is enough to compute Jf{m) for one m > A(/) — lct(/). 
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When it is known that the multipher ideal Jif^) is 0-dimensional, it is possible to 
bypass the elimination step (line [7] of Algorithm IS.lip in the following fashion. For a fixed 
monomial ordering > on -ft' [a;], we know that there are finitely many standard monomials 
(monomials not in the initial ideal m.yj{f'^)). Let h' G Q[s] be the polynomial produced 
by lines [8] and |9] of the above algorithm. A basis for the /^-linear relations amongst 
{a;"6'((j) I \a\ < d} modulo Jf{m) gives a basis Pd for the i^-space of polynomials in 
J^if^) up to degree d. By starting with d = and incrementing d until all monomials of 
degree d belong to m>{Pd-i), we obtain (P^) = J{f^) upon termination. 

Algorithm 5.12. Jif) = multiplierldealLA (/, c, cimax) 
Input: / = {/i, . . . , /,} c K[x], c E Q, d^ax e N. 

Output: the multiplier ideal J {f^) C A'[a;] , when it is generated in degrees at most d^ax- 

G,^{t^-f,\3 = l,...,r} U {d,^ + Y:'.=i'^dt, M = l,...,n}. 

m ^r- [max{c — lct(/), 1}] . 

if c — lct(/) is an integer and > 1 then 

m <(— m + 1 
end if 

G2 ^ starIdeal(Gi, w) U {/" | a G N'', \a\ = m} C Dy, where w assigns weight 1 to 
all df and to all dr ■ 

i? ^ I a Grohner basis of G2 w.r.t. any monomial order |. 
b ^ generalB(/, 1, m). 

b' ^ product of factors {s — c')"*-'^'-' ofb over all roots d ofb^^^ such that —d > c, where 

a{d) equals the multiplicity ot the root d . 

d < 1; P ^ C A'[a3] with > that respects degree. 

while P = ^ or (in> (P) does not contain all monomials of degree d and d < d^ax) 
do 

d ^ d + 1, 

{a\\a\< d, a;" ^ in>(P)}. 
Find a basis Q for the K-syzygies {qa)aeA such that 

J2q^mB{x''b{a)) = 0. 



P^PU{E„eA^aa;"|(g«)Gg}. 
end while 
return (P) . 

Notice that with dmax = 00 the algorithm terminates in case dim J'{f'') = 0. It also can 
be used to provide a A'-basis of the up-to-degree-dmax part of an ideal of any dimension. 
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6. Examples 

We have tested our implementation on the problems in [32]. In addition, this section 
provides examples from other sources with the theoretically known Bernstein-Sato poly- 
nomials, log-canonical thresholds, jumping numbers, and/or multiplier ideals; below is 
the output of our algorithms on several of them. 

The authors would like to thank Zach Teitler for suggesting interesting examples, some 
of which are beyond the reach of our current implementation. We also thank Takafumi 
Shibuta for sharing his script (written in risa/asir [23]), which is the only other existing 
software for computing multiplier ideals. 

A note on how to access Macaulay2 scripts generating examples, including the ones in 
this paper and some unsolved challenges, is posted at [3] along with other useful links. 

Example 6.1. When f = + y'^ + x^y"^, Saito observed that not all roots of bf{—s) are 
jumping coefficients [271 Example 4.10]. The roots of bf{—s) within the interval (0, 1] are 

9 11 13 7 17 9 19 

20' 20' 20' To' 20' To' 20' 
However, ^ is not a jumping coefficient of /. This can be seen in J/(l) from Theorem l5.10[ 
which has, among others, the primary components {s + ^, y, x) and (s -|- ^, |/, x). In fact. 



'C[x,y] 
x,y) 
x^,y) 



20' 

ifO<c<|, 
if A < c < 1^ 

20 — ^ ^ 20' 



if 13 < C < ^ 

20 — ^ 10 ' 

if X < c < 

10 — ^ 20 ' 

if IZ < c < ^ 

20 — ^ 1(1) 



(X 



xy,y^) 
xy,y'^) 



for all c > 1. 



20 
9 



10' 
19 



if JL < C < 
10 — ^ 20' 

if i < C < 1, 



and Jif) = if) ■ J{r') 

Example 6.2. We compute Bernstein-Sato polynomials to verify examples corresponding 
to [3ll Example 7.1]. The C[x, 2;]-ideal 

(/) = {x - z,y - z) n {^x - z,y - 2z) n (5?/ - x, z) 

defining three non-collinear points in has 

hf{s) = {s + \){s + 2f. 

In particular, its log canonical threshold is |. The multiplier ideals in this case are 



if < c < |, 
if I < c < 2, 



and J{f 



€[x,y,z\ 
Jx,y, z) 

if) ■ J^{f^~^) foi' all c > 2. On the other hand, the C[x, y, z]-ideal 
(g) = {y, z)n{x- 2z, y- z)n{2x- 3z, y- z) 
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defines three coUinear points in P^. Since 



the log canonical threshold of g is |. Here the multiplier ideals are 



C[x,y,z] 



if < c < |, 
if I < c < 2, 



and i7(g'^) = (g) ■ i7(g^ ''") for all c > 2. Thus, as Teitler points out, although g defines a 
more special set than /, it yields a less singular variety. 



Example 6.3. Consider / 



— — z'^){y'^ — z'^)z, the defining equation for a 



nongeneric hyperplane arrangement. Saito showed that | is a root of bf{—s) 
a jumping coefficient [261 5.5]. We verified this, obtaining the root 1 of bf{ 
multiplicity 3, as well the following roots of multiplicity 1 (including |): 

34256894 10 11 
7' 7' 3' 7' 7' 7' 7' S'Y'Y' 



but not 
-s) with 



Further, 



x,y,z) 

x,y,zf 

z,x) n {z,y)r\ 

y + z,x + z) r\ {y + z,x — z)r\ 
y — z,x + z)n{y — z,x — z) 
z,x) n {z,y)n 

y + z,x + z) n {y + z,x — z)n 
y — z,x + z) n {y — z,x — z)n 

■H 2 2 S S 2 2\ 

z^, yz , xz ,xyz, y", x , x y ) 



if < c < f , 
if I < c < ^ 



7 

if ^ 



7' 

- < C < - 



if I < C < f , 



if I < c < 1, 



and J{n = if) ■ J{r~^) for all c > 1. 



All examples in this section involve multiplier ideals of low dimension. In our experience. 
Algorithm 15.121 for a multiplier ideal of positive yet low dimension with a large value for 
o^max runs significantly faster than Algorithm 15.111 This is due to the avoidance of an 
expensive elimination step. 
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